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ABSTRACT 

In this paper, a robust technique is proposed for a single carrier frequency division multiple access system 
channel estimation in flat- fading non- Gaussian channels. Further, A new M-estimator is proposed for the 
robustification of the proposed detector. Simulation results are also provided to demonstrate the efficacy of the 
proposed detector over the least squares, Huber and Hampel based detectors in flat- fading non- Gaussian 
channels. 
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I. INTRODUCTION 
I. INTRODUCTION 

Wide demand on high data rates in wireless 
communication systems has arisen in order to 
support broadband services. The Third Generation 
Partnership Project (3GPP) Long Term Evolution 
(LTE) radio access standard provides peak data rates 
of 75 Mb/s on the uplink and 300 Mb/s on the 
downlink. In LTE standard orthogonal frequency- 
division multiple accesses (OFDMA) is used on the 
downlink. This supports different carrier bandwidths 
(1.25-20 MHz) in both frequency-division duplex 
(FDD) and time-division duplex (TDD) modes. In 
OFDMA each user is provided with a unique fraction 
of the system bandwidth. OFDMA combines 
scalability, multipath robustness, multiple-input 
multiple-output (MIMO) compatibility [1], thereby 
making it adaptive for wideband wireless 
accessibility. 

OFDMA, being sensitive to frequency offset and 
phase noise, accurate frequency and phase 
synchronization is needed. In addition, OFDMA is 
characterized by a high transmit PAPR, and for a 
given peak power limited amplifier this results in a 
lower mean transmit level. For these reasons. 



OFDMA is not well suited to the uplink 
transmission. Hence, LTE proposed, single carrier 
FDMA (SC-FDMA), also known as discrete Fourier 
Transform (DFT) pre coded OFDMA, for the uplink. 
PAPR reduction in SCFDMA is motivated by a desire 
to increase the mean transmit power, improve the 
power amplifier efficiency, increase the data rate, 
and reduce the bit error rate (BER) and energy 
consumption [4]. SC-FDMA ensures high data rate 
transmission, utilizing single carrier m frequency 
domain equalization. In algorithm is proposed for 
LTE uplink channel impulse response (CIR) 
knowledge of the channel. 

II. CHANNEL MODEL 

At the transmitter side, a baseband modulator 
transmits the binary input to a multilevel sequences 
of complex numbers ml(q) in one of several possible 
modulation formats including binary phase shift 
keying (BPSK), quaternary PSK (QPSK), 8 level PSK 
(8-PSK), 16-QAM, and 64-QAM. These modulated 
symbols perform an N-point discrete Fourier 
transform (DFT) to produce a frequency domain 
representation. The DFT followed by IDFT in a 
distribution-FDMA (DFDMA) or localization-FDMA 
(LFDMA) subcarrier mapping setup operates as 
efficient implementation to an interpolation filter. In 
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distributed subcarrier mode, the outputs are 
allocated equally spaced subcarriers, with zeros 
occupying the unused subcarrier in between. While 
in localized subcarrier mode, the outputs are 
confined to a continuous spectrum of subcarriers. 
Except the above two modes, interleaved subcarrier 
mapping mode of FDMA (IFDMA) is another special 
subcarrier mapping mode [12], [13]. The difference 
between DFDMA and IFDMA is that the outputs of 
IFDMA are allocated over the entire bandwidth, 
whereas the DFDMA' s outputs are allocated every 
several subcarriers. The response of the wide-sense 
stationary uncorrelated scattering (WSSUS) fading 
channel can be represented as 



/-I 



j=0 (1) 

where fading channel coefficients Wi(t) are the wide 
sense stationary i.e. W](t) = w(m; j), uncorrelated 
complex Gaussian random paths gains at time 

instant t with their respective delays T j, where w(m; 

j) is the sample spaced channel response of the 1*^ 
path during the time k, and S(^) is the Dirac delta 

function. Based on the WSSUS assumption, the 
fading channel coefficients in different delay taps are 
statistically independent. And has an autocorrelation 
function given by 



E {^(k, j)w^ (n, j) = ( j) Jo J;/^7} (k-n)_ 



(2) 



Where w (n; j) is a response of the 1*^ propagation 
path measured at time n, cr^(/) denotes the power 

of the channel coefficients, fd is the Doppler 
frequency in Hz, Tf is the symbol duration in 

seconds, and Jq (•) is the zero order Bessel function 
of the first kind (4). 

III. PROPOSED ALGORITHM 

The signal s (k) is transmitted via a time-varying 
channel w (k), and corrupted by an observation noise 
n (m) before being detected in a receiver. The signal 
received at time index k is 



r{k) = s^{k- 1) (m) + + - /) (k) + n{k) 

I 

7=1 

= S'{k)W{k) + n{k) 

Where s j(k — j) , j = 1,2, . . . , / are transmitted 

signal vectors at time m, 1 is the distinct paths from 
transmitter to the receiver, w(m) is the channel 
coefficients at time m, and n(k) is the noise with zero 

mean and variance . After processing some 

intermediate steps (synchronization, remove CP, 
DFT, and de mapping), the decision block 
reconstructs the detected signal to an approximate 
modulated signal and its phase. The output y(k) of 
the adaptive filter is expressed as 

y(k) = d^(k- 1)/Zj (k)-\- -{-d^ik- l)hi (m) 

= j^d.{k-j)h.{k) 

7=1 

= D^{k)h{k) 

Where d j (k — j), j = 1,2, . . . , / are detected signal 
vectors at time 

k,.. D(k) = diag |?j (/: - 1) , ^2 ~ 2) , . . . , (/: - /) 

In this problem formulation, the ideal adaptation 
procedure would adjust W](k) such that 
W](k)=h](k)as ^ ^ oo . In practice, the adaptive filter 
can only adjust w(k) such that y(k) closely 
approximates desired signal over time. There- fore, 
the instantaneous estimated error signal needed to 
up- date the weights of the adaptive filter is 

j(k) = p{k)e^(k)e{k) 
e(k) = r{k) - y(k) 

= r{k) - (k)h{k) 

This priori error signal, e(k) is used to minimize the 
estimator error by adaptive updating of filter 

weights. w(0 = [Wq(0,.— Wa^_i(0] is the channel 

noise vector during the i^^ symbol interval. w(i) are 
assumed to be independent and identically 
distributed random variables with non- Gaussian 
distribution. [ ] 
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(l-s)/(0,af) + s/(0,a^) 



(12) 



Several different classes of robust methods (such as 
M, L, R estimators, and the least median of squares 
method) exist [7-9]. In this paper, M-estimator based 
robust clustering technique is proposed. The M- 
estimators try to reduce the effect of outliers by 
replacing the squared residuals by less rapidly 

increasing function p of the residuals, yielding 



nun 

c 



N 



(3) 



Where p is a symmetric, positive-definite function 

with a unique minimum at zero, and is chosen to is 
less increasing than squared function? Let 

C = [Cp <^yt]be the parameter vector to be 

estimated. The M-estimator of c based on the 

function /^(^ ) is the vector C which is the solution 
of the following k equations: 



dr 



(4) 



dCj 



dp(x) 

where the derivative (//{x) = is called the 

dx 



influence function. The influence function 
l//(x) measures the influence of a datum on the value 

of the parameter estimate. In this paper, a new M- 
estimator is proposed for robustifying the k-means 
clustering algorithm. The penalty function and the 
influence functions of the proposed M-estimator are 
given by [10] (also see Fig 1). 
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where d is any constant. 
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The choice of the constants a = kv^ and b 

= 2kv depends on the robustness measures 
derived from the influence function. 




Fig. 1(a) 




Fig.l (b) 

Fig. 1 (a) Penalty functions and (b) influence 
functions of the Proposed M-estimator. 

IV. SIMULATION RESULTS 

The performance of the proposed algorithm is 
compared with the Least Squares, Huber and 
Hampel based detectors in a Rayleigh fading 
environment and corrupted with non-Gaussian 
noise. The simulation parameters are listed in tablel. 
The BER is a significant performance parameter for 
quality measurement recovered data in wireless 
communication effect of the proposed algorithm in 
terms of BER compared with existing estimators. It is 
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evident from these results that the proposed 
algorithm outperforms the Least Squares, Huber and 
Hampel based detectors in non-gaussian 
environment. 
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Fig. 2 Probability of error versus SNR for the 
considered detectors in Gaussian noise. 
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Fig. 3 Probability of error versus SNR for the 
considered detectors in non-Gaussian noise. 
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Fig. 4 Probability of error versus SNR for the 
considered detectors in non-Gaussian noise. 



V. CONCLUSIONS 

In this paper, a robust detection technique is 
proposed to estimate channel parameters when the 
signal is corrupted with non-Gaussian noise. A new 
M-estimator is proposed. Simulation results are also 
provided. These simulation results prove that the 
proposed detector outperforms Least Squares, Huber 
and Hampel based detectors in non-Gaussian flat- 
Fading channels. 
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